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Abstract 

Taking into account the gluon shadowing and the J I if/ nuclear absorption, we evaluate the Cold Nuclear Matter 
(CNM) effects on J /if/ production in pPb and PbPb collisions at the design LHC energy for PbPb collisions, s/snn = 
5.5 TeV. In view of the good agreement between the yield prediction using the LO CSM and the ALICE data at 7 TeV, 
we employ the latter model to use a complete description of the kinematics of the underlying 2 — > 2 partonic process, 
namely g + g — > J/if/ + g. We observe a large J/ if/ suppression in /?Pb and PbPb collisions due to the strong shadowing 
of the gluon distribution in the lead ion at the corresponding gluon momentum fractions. This suppression from CNM 
effects has a strong rapidity dependence which may compensate, partially or completely, that of the predicted charm 
recombination in PbPb collisions. 
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1. Introduction 

Quark Gluon Plasma (QGP), a deconfined state of QCD matter, is expected to be produced by Relativistic nucleus- 
nucleus (AB) collisions at high enough densities or temperatures. The J/if/ meson should be [1] sensitive to Hot 
and Dense Matter (HDM) effects, through processes like the colour Debye screening of the cc pair. A significant 
suppression of the J /if/ yield was observed at SPS energy by the NA50 experiment [2], and at RHIC by the PHENIX 
experiment in CuCu [3] and AuAu [4] collisions at y/s^N = 200 GeV. In 2010, data have been taken at the LHC in 
PbPb collisions at y^Aw = 2.76 TeV. After the long LHC shutdown planned in 2012, LHC will hopefully run at its 
design energy with PbPb collisions at -\Js NN = 5.5 TeV. These results will give the possibility to test the available 
models at a new energy scale. However, the interpretation of the results obtained in AB collisions relies on a good 
understanding and a proper subtraction of the Cold Nuclear Matter (CNM) effects which are known to already impact 
the J/if/ production in proton (deuteron)-nucleus (pA or dA) collisions, where the deconfinement cannot be reached. 
Two CNM effects are of particular importance [5]: (i) the shadowing of the initial parton distributions (PDFs) due to 
the nuclear environment, and (ii) the breakup of cc pairs after multiple scatterings with the remnants of the incident 
nuclei, referred to as the nuclear absorption. In our previous works [6-9], we developed an exhaustive study of these 
effects at yjs^N - 200 GeV and we confronted our results to the measurements from PHENIX [10]. It is our purpose 
here to extend our results to /?Pb and PbPb collisions at the energy yjs^N - 5.5 TeV. 

We have shown in earlier studies [6-9] that J /if/ partonic production mechanism affects both the way to compute 
the nuclear shadowing and its expected impact on the J/ iff production. Most studies on the J /if/ production in hadronic 
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collisions assume that the cc pair is produced by a 2 — > 1 partonic process where both initial particles are two gluons 
carrying some intrinsic transverse momentum kj. The sum of the gluon intrinsic kj is transferred to the cc pair, 
thus to the J/if/ since the soft hadronisation process does not significantly alter the kinematics. This is supported by 
the picture of the Colour Evaporation Model (CEM) at LO (see [11] and references therein) or of the Colour-Octet 
(CO) mechanism at a] [12]. Thus, in such approaches, the transverse momentum Pj of the J /if/ comes entirely 
from the intrinsic kj of the initial gluons. This is not sufficient to describe the Pj spectrum of quarkonia in hadron 
collisions 1 [11]. As a consequence, J /if/ production at low and mid Pj likely proceeds via a 2 — ♦ 2 process, such as 
g + g — > J/i// + g, instead of a 2 — > 1 process 2 . This amounts to the bulk of the J /if/ production cross section. For 
instance, as illustrated by Fig la, the yield predicted by the LO CSM [18] reproduces correctly the PHENIX, CDF 
and ALICE measurements without resorting to any colour-octet mechanism nor parameter twisting. Consequently, 
one is entitled to consider that the former 2 — > 2 kinematics is the most appropriate to derive CNM effects at RHIC, 
and to provide predictions at LHC energy. In this work, we shall focus on the CNM effects expected at the energy 



ysjvjv = 5.5 TeV in the extrinsic scheme. The article is organized as follows: in section 2, we will describe our model 
and in section 3, we will present and discuss our results. 

2. Our approach 

To describe the J /iff production in nuclear collisions, our Monte Carlo framework [6, 19] is based on the proba- 
bilistic Glauber model. The nucleon-nucleon inelastic cross section at y/s NN = 5.5 TeV is taken to be cr NN = 
70 mb [20] and the maximum nucleon density to be po = 0.17 nucleons/fm 3 . We also need to implement the partonic 
process for the cc production model that allows us to describe the pp data and the CNM effects. 

For Pj >2 - 3 GeV, most of the transverse momentum of the quarkonia should have an extrinsic origin, i.e. the 
J/iff's Pj would be balanced by the emission of a recoiling particle - a hard gluon - in the final state. The J /if/ would 
then be produced by gluon fusion in a 2 — > 2 process. This emission, which is anyhow mandatory to conserve C- 
parity, has a definite influence on the kinematics of the J/ifi production. Indeed, for a given J/ifi momentum (thus for 
fixed rapidity y and Pj), the processes discussed above, i.e. the intrinsic g + g — > cc — > J/if/(+X) and the extrinsic 
g + g — > J/if/ + g, will proceed from initial gluons with different Bjorken-x on the average. Therefore, they will be 
affected by different shadowing corrections. 

In the intrinsic scheme, the measurement of the J /if/ momentum in pp collisions completely fixes the longitudinal 

momentum fraction of the initial partons: jq 2 = d|-= ex P (±v) = x® 2 (y, P T ) with m T - -^M 2 + P T , M being the J /if/ 
mass. On the contrary, in the extrinsic scheme, the knowledge of the y and Pj spectra is not sufficient to determine x\ 
and X2- Actually, the presence of a final-state gluon introduces further degrees of freedom, allowing several (x\ , X2) for 
a given set {y, P T ). The four-momentum conservation results in a complex expression of xi as a function of (x\ , y, P T ): 



kinematics determines the physical phase space, models are anyhow mandatory to compute the proper weighting of 
each kinematically allowed (x\, X2). This weight is simply the differential cross section at the partonic level times the 
gluon PDFs, i.e. g{x\,pj )g(x2,fiF ) dcr gg ^ji^ + ^/dy dP T dx\dxi- In the present implementation of our code, we are able 
to use the partonic differential cross section computed from any theoretical approach. In this work, we shall use the 
Colour-Singlet Model (CSM) at LO at LHC energy, shown to be compatible (see Fig. la) [14, 18] with the magnitude 
of the P7- -integrated cross-section as given by the PHENIX pp data [21], the CDF pp data [22] and the LHC pp data 
at = 7 TeV. 



To obtain the J/ifi yield in pA and AA collisions, a shadowing-correction factor has to be applied to the J/if/ yield 
obtained from the simple superposition of the equivalent number of pp collisions. This shadowing factor can be 



'in addition, recent theoretical works incorporating QCD corrections or .v-channel cut contributions have emphasized [13-15] that the Colour- 
Singlet (CS) mediated contributions are sufficient to describe the experimental data for hadroproduction of both charmonium and bottomonium 
systems without the need of CO contributions. Furthermore, recent works [16] focusing on production at e + £>~ colliders have posed stringent 
constraints on the size of CO contributions, which are the precise ones supporting a 2 — > 1 hadroproduction mechanism [11]. 

- One may also go further and consider more than two particles in the final state, as expected from the real-emission contributions at NLO and 
NNLO [13]. It is clear from the yield polarisation [17] that these contributions start to dominate for Pj above 1 - 2m c . The effect of more partons 
in the final state is to increase the difference between the results obtained in both schemes. However the implementation of NLO and NNLO codes 
in a Glauber model with an inhomogeneous shadowing is not yet available. 



X2 = 




. Equivalently, a similar expression can be written for x\ as a function of (x2,y, Pt)- Even if the 
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expressed in terms of the ratios Rf of the nuclear Parton Distribution Functions (nPDF) in a nucleon belonging to a 

nucleus A to the PDF in the free nucleon: Rf(x, Q ) = A ^ tldm(x Ql) , i = q,q,g . The numerical parameterisation 

of Rf (x, Q 2 ) is given for all parton flavours. Here, we restrict our study to gluons since, at high energy, the J /if/ is 
essentially produced through gluon fusion [11]. Several shadowing parametrisations are available [23, 24]. In the 
following, we shall restrict ourselves to EKS98 [23], which is very close to the mean in the current evaluation of 
the uncertainty [24] on the gluon nPDF and exhibits a moderate antishadowing. We postpone the propagation of the 
uncertainty on the gluon nPDF to the CNM effects evaluated at LHC energy for future studies. 

The second CNM effect that we are going to take into account concerns the nuclear absorption. In the framework 
of the probabilistic Glauber model, this effect is usually parametrised by introducing an effective absorption cross 
section <x abs . It reflects the break-up of correlated cc pairs due to inelastic scattering with the remaining nucleons from 
the incident cold nuclei. The value of cr abs is unknown at LHC. At high energy, the heavy state in the projectile should 
undergo a coherent scattering off the nucleons of the target nucleus [25], in contrast with the incoherent, longitudinally 
ordered scattering that takes place at low energies. As argued in [26, 27], this should lead to a decrease of <x abs with 
increasing ^[snn- The systematic study of many experimental data indicate that cr abs appears either constant [28] or 
decreasing [29] with energy. Hence, we can consider our estimates [8, 9] of cr abs at RHIC energy as upper bounds for 
the value of cr abs at LHC. We choose three values of <x abs that should span that interval (cr abs = 0, 1.5,2.8 mb). 




Figure 1: (a) dcr^ ct / dy\ Y= o x Br from gg fusion in pp collisions for yfs from 200 GeV up to 14 TeV compared 
to the PHENIX [21], CDF [22] and ALICE data, (b)-(c) J /if/ nuclear modification factor versus y in pPb and PbPb 
collisions at = 5.5 TeV, using EKS98 [23] gluon shadowing parametrisation and three values of cr a i, s (from 

top to bottom: 0, 1.5,2.8 mb) in the extrinsic scheme. For PbPb collisions, the y-dependence is shown for various 
centrality selections. 

In the following, we present our results for the J /if/ nuclear modification factor due to CNM effects in the extrinsic 
sheme in pPb and PbPb collisions at V^aw = 5.5 TeV: R AB = dN J J* /{N ma )dN J J^, where dN J /% (dN J J*) is the 
observed J/if/ yield in AB = pPb, PbPb (pp) collisions and (N co n) is the average number of nucleon-nucleon collisions 
occurring in one pPb or PbPb collision. Without nuclear effects, Rab should equal unity. 

In Fig. lb, we show R p pb versus y. The curve with no absorption allows us to highlight the strong rapidity 
dependence of the shadowing. We can also notice that the shadowing alone should already be responsible for a quite 
large amount of J /if/ suppression, up to 36 % at y = 6. This is expected due to the very small x-region in the gluon 
nPDF that becomes accessible at LHC energy (down to 10~ 5 ). At backward rapidity, we are in the antishadowing 
region, with the antishadowing peak at y ^ -5. Fig. lc shows that the y-dependence of Rp^pb for the centrality bins 
used by the ATLAS Collaboration [30] is similar for all the bins, with a dip at mid-y. This shape is the opposite of 
the one obtained at RHIC energy [6, 8], with a peak at mid-y. One can also notice that Rpbpb decreases at very large 
rapidity showing that one has gone beyond the antishadowing peak. Here, /?p b p b is systematically smaller at mid->' 
than at forward-y. This is also illustrated on Fig. 2, with the centrality dependence of /?p b p b for two regions in y. 
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Figure 2: (Color online) J/ip nuclear modification factor, ^pbpb, in PbPb collisions at ysjyjv = 5.5 TeV versus Np^, 
using EKS98 [23] gluon shadowing parametrisation and three values of the nuclear absorption cross section in the 
extrinsic scheme. Rp\,pb is shown for two different experimental acceptances in rapidity, |y| < 0.9 and \y\ e [2.5,4]. 



This behaviour of the CNM effects may partially - or completely - compensate the opposite effect expected from cc 
recombination, with a maximum enhancement at y — 0. Overall, one may observe a Rp^pb rather independent of y 
resulting of two y-dependent effects. 
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